Decay widths of charmonia in a hot equilibrated medium 
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We investigate the properties of charmonia in a thermal medium, show- 
ing that with increasing temperature the decay widths of these mesons be- 
have in a non-trivial way. Employing a potential model with interaction po- 
tential extracted from thermal lattice QCD calculations of the free-energy 
of a static quark- antiquark pair, we study some decay processes in the 
crossover region. We find that at temperatures T ~ T c the decay widths 
of the J/ty that depend on the value of the wave function at the origin 
are enhanced with respect to the values in vacuum. In the same tempera- 
ture range the decay width of the process \c.j — * J + 7 is enhanced by 
approximately a factor 6 with respect to the value in vacuum. At higher 
temperatures the charmonia states dissociate and the widths of both decay 
processes become vanishing small. 

PACS numbers: 12.38.-t, 12.38.Aw, 12.39.Pn, 13.20.Gd, 25.75.Nq 

1. Introduction 

The heavy-ion collision program aims to identify and characterize the 
properties of deconfined quark matter at high temperatures. Many valuable 
probes of the properties of the medium produced in a heavy-ion collision are 
available, which include jets, electromagnetic signals and heavy quarkonia 
states (QQ) [TJ. In particular much work has been devoted to understand 
how the presence of deconfined quarks and gluons may affect the binding of 
quarkonia [21 El H] . 
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The analyses of the binding energies of mesonic states as a function 
of the temperature have been carried out by various authors using poten- 
tial models, with interaction potentials extracted from lattice QCD cal- 
culations H |j. These analyses show that quarkonia dissociate at 
temperatures close to the critical temperature of QCD, T c . Analogous re- 
sults have been obtained by studying the correlation functions of charmonia 
above deconfinement [9]. 



In the present study [10] we are interested in determining how the decay 

widths for the leptonic, hadronic and radiative channels are influenced by 
the presence of the thermal medium. In order to determine the decay widths 
we first evaluate the radial wave function of the pertinent charmonia states 
employing the non-relativistic Scrhodinger equation (and for comparison 
the relativistic Salpeter equation) with interaction potential extracted from 
thermal lattice QCD simulations. 

The remarkable result of our analysis is that some decay widths change 
drastically close to T c . In the transition region we find that the decay widths 
of the J/ty that depend on the value of the wave function at the origin change 
by approximately a factor 2 with respect to the corresponding values in 
vacuum. We also investigate the properties of the Xc meson and in particular 
of the Xcj —> J/^ + 7 transition. This radiative decay contributes to the 
total inclusive J/^> production by a fraction of about 0.3, as determined 
in proton proton and n + tt~ data [11]. We find that this J/^> production 
mechanism is enhanced close to the critical temperature by approximately 
a factor 6 with respect to the value in vacuum. 



Potential models have been quite successful in describing the properties 
of heavy quark bound states in vacuum (see e.g. [E]). Our key assumption 
is that at any temperature T, the interaction between heavy quarks can be 
approximated by an instantaneous potential, V(r,T), where r is the radial 
coordinate. 

We shall study bound states of heavy quarks by using the non relativistic 
Schrodinger equation 



where Mq is the constituent mass of the heavy quark; ipi and E{ represent 
the wave function and the mass of the corresponding QQ state, respectively. 
Comparing the results of this analysis with those obtained using the Salpeter 
equation for the S— wave states 




2. The method 
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we find excellent agreement between the outcomes of the two methods, 
meaning that for these mesonic states relativistic corrections are small. 

There are some controversies about how to extract the potential V(r, T) 
from the color-singlet free-energy F\(r,T) (see [BlElllS])- I n principle, the 
internal energy U± is obtained subtracting the entropy contribution from 
the free-energy 

^F^ 

U 1 = F 1 -T^. (3) 

However, this subtraction procedure has been put to question, one of 
the reasons being that around the critical temperature, the potential U\ is 
more attractive than the potential at zero temperature. One possible inter- 
pretation of this result is that at non-vanishing temperatures there can be 
additional interactions between the two static quarks [13] . Below T c this 
effect should be due to hadrons and could be related to the string "flip-flop" 
interaction [II]; above T c it should be related to the antiscreening properties 
of QCD. In the following and in agreement with Ref. [13] we shall use the 
internal energy in Eq. ([3]) as the interaction potential and show the results 
obtained using the free-energy only for comparison. As we shall see, us- 
ing Ui, we obtain dissociation temperatures for charmonia consistent with 
lattice results obtained with the Maximum Entropy Method [9]. Employ- 
ing the free-energy as interaction potential gives much smaller dissociation 
temperatures. 

We shall consider the free-energy proposed in [15] employing the formal- 
ism of the Debye-Hiickel theory, with the parameterization of Ref. [16] , 

OL 

(exp(— x)+x) , 

r 

where is the modified Bessel function, x = fir, while the functions 

fi = fi(T) and n = k(T) are determined by fitting the lattice data. Once 
these two functions are fixed, this parameterization of the free-energy is 
in excellent agreement with lattice data for T < 2T C . At short distances 
the static quark-antiquark free-energy is normalized in such a way that it 
reproduces the free-energy at zero temperature 

01 

F(r,T = 0) = +ar, (5) 

r 

where a is the string tension and a is the Coulomb coupling constant. The 
values of these two parameters can be phenomenologically fixed; we shall 
take a = 0.16 GeV 2 , a = 0.2 and with m c = 1.28 GeV one obtains Mj/^ ~ 
3.11 GeV and M v ~ 3.43 GeV. 



Fi(r,T) 



a 
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3. Masses and Decay widths 

In order to determine the decay widths it is necessary to compute the 
wave function and the mass of the charmonia states. For I = states, 
they are determined solving both the Schrodinger equation and the Salpeter 
equation through the Multhopp method while for P— wave states we only 
solve the Schrodinger equation. 

We have considered various decay processes. Annihilation processes of 
charmonia can be viewed as two-stage factorized processes. The first process 
consists in cc annihilation into gauge bosons. The annihilation takes place 
at the characteristic distance r of order l/m c , so for a non-relativistic pair 
r — > and the annihilation amplitude is proportional to the wave function at 
the origin. Then, the produced gauge boson decays or fragments into leptons 
and hadrons. In vacuum it is assumed that the inclusive probability of the 
latter process is equal to one. At finite temperature we assume that such a 
factorization persists. Therefore, the width of hadronic and leptonic decays 
will depend on temperature exclusively through the radial wave function at 
the origin and the mass of the meson [12] and can be expressed as 

T A oc Y l+l - oc T( 3 Si -> 3g) oc '^^f . (6) 

We also consider the radiative transitions XcJ — ¥ J/\P + 7 from the 3 Pj 
levels to the 3 Si state with rate 

T B oc T XcW/ * +7 oc (2J + 1)|/ PS | 2 , (7) 

where Ips is the overlap integral between the radial wave functions of the 
corresponding states. 



3.1. Decay widths at different temperatures 

Now we show the results for dissociation temperatures and decay widths 
obtained using as QQ potential the internal energy and for comparison the 
free-energy. 

Using the free-energy as interaction potential, the J/ip dissociates at 
temperatures close to 1.2 T c , while the Xc dissociate at about 0.95 T c . Using 
the internal energy as potential we find that the J/ ip dissociates at approx- 
imately 1.8 T c , while the \ c dissociates at about 1.15 T c . Using the internal 
energy one obtains higher dissociation temperatures in agreement with the 
fact that U\ is more attractive than F\. 

In Fig. Q] we report the values of the decay width for the process in 
Eq. ([6]), left panel, and for the process in Eq. (JJJ), right panel, as a function 
of the temperature. Employing F\ one finds a monotonic decrease of the 
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width that can be viewed as due to a decrease of the screening length of 
the potential. On the other hand, when using Ui we find that across the 
transition region both decay widths become large. For the decays in Eq. ([6]) 
the width is about a factor 2 larger than in vacuum. For the radiative 
decay in Eq. ([7]), the ratio between the width in the thermal medium and 
in vacuum reaches a factor 6 across the transition region. This effect is 
not due to the screening of the potential, but is rather connected with the 
"cloud-cloud" interaction [13J. 

In order to have a rough estimate of the effect, let us suppose that the 
Xc is produced in the early stage of the heavy-ion collision and then travels 
for about 4 fm in the thermally equilibrated medium. If the temperature 
is sufficiently low, so that the decay width is approximately the same as in 
vacuum, one has that less than 1% of Xc decay. On the other hand, for 
temperatures close to T c , one should observe that about 5% of Xc decay. 
Since for temperatures close to T c one can neglect interactions of charmonia 
with in-medium hadrons and gluons [2] the radiative decay process should 
give a sizable contribution to the total decay width of the Xc- As regard 
the inclusive width of the J/^, one should consider that at temperatures 
below T c , the process in Eq. © gives approximately the inclusive width. 
However, at larger temperatures the dominant contribution should be due 
to the interaction with in-medium partons [2]. 




T/T T/T 

c t 



Fig. 1 . (color online) Charmonia decay widths given in Eq. ([6|) (left panel) and in 
Eq. ([7]) (right panel). We have used iq, full blue lines, and Lq, dashed red lines, as 
potentials with the paramterization of Eq. (j4|) , and normalized the widths to the 
vacuum value. 

The largest uncertainty in our calculations resides in the extraction of 
the potential from lattice QCD calculations and pertinent parametrizations 
to numerically evaluate the entropy. In order to test the robustness of our 
results one can employ a different parametrizations of the internal energy. 
Using the expression of U\ (r, T) reported in Ref . [6] we find approximately 
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the same results reported here. 

The analysis reported here can easily be extended to bottomia, moreover 
it would be of some interest to include D-mesons. 
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